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We present an efficient method for subrecoil cooling of neutral atoms by applying Raman cooling 
in 2D to a four-level tripod-system. The atoms can be cooled simultaneously in two directions using 
only three laser beams. We describe the cooling process with a simple model showing that the 
momentum distribution can be rapidly narrowed to velocity spread down to O.lt'rec, corresponding 
to effective temperature equal to O.OlTrec- This method opens new possibilities for cooling of neutral 
atoms. 

PACS numbers: 32.80.Pj 



I. INTRODUCTION 

Cold atoms are a multipurpose tool for both fundamen- 
tal observations and technological innovations. Cooling 
and trapping of neutral atoms has led to many discov- 
eries such as Bose-Einstein condensation Jil, atom in- 
terferometry [2] , atomic nanofabrication |3|, 1^ and im- 
proved atomic clocks [E*, ^ . Evaporative cooling in mag- 
netic traps is an efficient tool for reaching high densities 
at low temperatures , necessary for the observation of 
quantum many-body effects in dilute gases. It has been 
applied with success also to optically trapped gases [H, 
as well as to the cooling of fermionic atoms via sympa- 
thetic cooling [9]. However, for many applications it is 
not necessary or even desirable to reach quantum degen- 
eracy. The external fields and atomic interactions can be 
harmful for ultraprecise measurements and atomic clocks, 
for which low temperatures coupled to low densities but 
large numbers of atoms are needed for a good signal-to- 
noise ratio and unperturbed atomic transition frequen- 
cies. Evaporative cooling is rather wasteful on atoms, 
and it requires also a large cross-section of elastic colli- 
sions and strongly confining traps with velocity-selective 
output coupling. For quantum information purposes one 
has to target single atoms, and then cooling without col- 
lisional thermalization is needed. Finally, the current 
interest in nanomechanics has stirred renewed interest 
on laser cooling |Tl| . We present an efficient Raman 
cooling method, which allows one to reach subrecoil tem- 
peratures rapidly with purely optical means in 2D, with 
a simple pulse setup. The method is applicable to all 
densities and even single atoms, and can be combined 
with sympathetic cooling if necessary. 

The original Raman cooling idea was presented and 
verified experimentally by Kasevich and Chu in 1992 at 
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Stanford for Na atoms in ID [l2|. The setup consists of 
a three-level A-type atom. A Raman pulse (two contra- 
propagating beams) moves an atom from one state to 
another and at the same time gives it a velocity change 
of 2'L'rec = 2hk/m. The pulse duration and detunings are 
selected so that the two-photon Raman process is sup- 
pressed for atoms with velocity near ^' = 0. The central 
frequency and duration of the subsequent pulses is ad- 
justed to target different velocity groups while still avoid- 
ing V = {). As atoms are optically pumped back to the 
original state after each pulse, they eventually accumu- 
late at velocities near v = {). The original work was later 
extended to 2D and 3D cooling [13], without reaching 
subrecoil temperatures. The ID scheme has been used 
in attempts to reach high 3D phase-space densities in 
optical traps [13, IHI . The use of a dipole trap allowed 
subrecoil 3D cooling by ergodicity of the motion in the 
inverted pyramidal trap fl3 or the simple anisotropy of 
a 3D harmonic trap Il7j while applying the Raman 
cooling only in ID. Similar work with Cs atoms has been 
done at ENS, Paris Whereas the original scheme 

used Blackman pulses of finite spectral width, it was ac- 
tually shown with Levy flight simulations that for 2D and 
3D the square pulses are quite efficient [18]. 

The key issue is that in the original scheme the di- 
rect Raman cooling in 2D requires four Raman beam 
pairs and can be quite cumbersome. The 2D scheme with 
square pulses was applied to Cs atoms at NIST, Gaithers- 
burg, giving temperatures down to O.lSTrec [19], where 
Tree is the atomic recoil temperature. This is actually 
the only demonstrated case of reaching subrecoil tem- 
peratures with true 2D Raman cooling. We show that in 
the tripod system, one can go down in theory by an or- 
der of magnitude or even more, which means going down 
to Y^Trec- The scheme can also be used in transversal 
cooling of atomic beams, for example Ne* [2Q|, [21|. 

The structure of this presentation is such that we first 
derive the basic equations for the velocity changes in- 
duced by one cycle of Raman pulses in the 2D tripod 
configuration. We show that in the limit of small ef- 
fective interaction between ground states one can obtain 
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The interaction of a four-level tripod-type atom with 
the laser beams is described by the Schrodinger equation 

The atomic Hamiltonian and the probability function in 
the rotating-wave approximation (RWA) and in the res- 
onant approximation are given by 

/-Ai Ge-'^^\ 

-A2 Ge-'^y 
-A3 Ge'^' 
\Ge'^^ Ge'^y Ge''^^ Q 

(2) 

where Q is the detuning of driving fields from the upper 
state, A^ = ujrn — ^4m — ^ {m = 1,2,3) are detunings 
from ground states |m) (see Fig. [Hb)). Without losing 
generality we have assumed that the Rabi frequencies 
Gj = \d4,jEj/h\ {j = 1,2,3) equal a real magnitude G. 

The probability amplitudes in the momentum-space 
approach are obtained by using the Fourier transform 
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FIG. 1. (a) An atomic ensemble interacts with three run- 
ning optical waves, (b) The level diagram of a four-level 
tripod-type atom. The cr+ and cr_ circularly polarized contra- 
propagating waves couple transitions |1) ^ |4) and |3) o |4), 
respectively, and the 7r-polarized wave couples the |2) ^ |4) 
transition. 



simple solutions for the equations. Then we extend the 
treatment to arbitrary interaction strength. Finally, we 
demonstrate the efficiency of the method by solving the 
model numerically for an appropriate initial momentum 
distribution, including both the Raman cycles as well as 
the optical pumping between them. 



II. 2D RAMAN COOLING IN A TRIPOD 
CONFIGURATION 

A. Basic equations 

Consider an atomic ensemble in the region of interac- 
tion with laser configuration that consists of three run- 
ning optical waves shown in Fig. [IJa). Two contra- 
propagating cr+ and cr_ circularly polarized waves are 
directed along the z direction, a 7r-polarized optical wave 
propagates along the y direction. The positive frequency 
part of the electric field can be written as 
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where k is the wave number same for each optical wave. 
The first term corresponds to a cr+ circularly polarized 
wave of frequency uji and the third one corresponds to a 
(J- circularly polarized wave of frequency uj^. The second 
term describes a 7r-polarized wave of frequency 002. 
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where p^, Py are projections of the atomic momentum on 
the axes Oz, Oy in units of hk. For variables 



hi = ai{pz - ^,Py), h = a2{pz,Py - 1), 

^3 = asiPz + l,Py), h = CLA{Pz,Py), 

the Schrodinger equation ([1]) can be written as 

ihi = {{pz - if + Si) hi - gh4, 

ih2 = {pI + {py - If + S2) h2 - gh^, 

ibs = {{Pz + 1)' + ^3) ^3 - gb4. 

ib^ = {pI + - A) 64 - gbi - gb2 - gbs, 

where the dimensionless time, detunings and the Rabi 
frequency are given by 



(3) 

(4a) 
(4b) 
(4c) 
(4d) 



ojr' ujr' 



(5) 



and ujji = hk^ /2M is the recoil frequency. 

The upper state |4) can be adiabatically eliminated in 
the case of Raman transitions between the ground states 
when A ^ g. In this case, we assume that 64 ~ and 
obtain from Eq. (|4d|) expression 
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Terms p^, Py were dropped here as values negligible in 
comparison with A. The substitution of Eq. (|6]) into 
Eqs. (j4|) leads to equations 
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where a = ^^/A is the Rabi frequency of the two-photon 
resonance. 

If the atomic population is initiahy concentrated in 
state 1 2) with corresponding amphtude it fohows from 
Eqs. (O that during the interaction 
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To solve Eqs. ([7]) one can consider some special cases. 
These are useful for gaining insight about the parameter 
choices such as pulse duration that optimize the efficiency 
of the approach. 



Short-time interaction 
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Let us consider the case of the short-time interaction 
when ar <C 1. In such a case, variables 61,63 ^ and 
therefore 61,63 vanish from the right-hand side of Eqs. 
([7]). After that, Eq. (|7b|) can be solved separately, and 
the substitution of 62 into Eqs. ([7a|) and ([7c|) leads to the 
solutions 
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where the denominators are 



S2)-Pz -^Vy, D2 = -{S3 



By substituting Eqs. ([9]) into Eq. ([8]), we get 
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Then the atomic populations in ground states, using 
Eqs. (j3j), ^ and (pTj) . can be written as 
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Let us consider atoms originally prepared in state |2) 
after they interact with a pulse of Raman beams. Cor- 
responding approximate forms of the populations in the 
momentum-space approach are given by Eqs. (p!2]) - (p!4|) . 
The population transferred to state |1) through Raman 
transition |2) ^ |1) is shown in Fig. Efa), while the pop- 
ulation transferred to state |3) through Raman transition 
1 2) ^ 1 3) is shown in Fig. Efc). The removed popula- 
tion of state 1 2) represents a cross as shown in Fig. Efb). 
The width of the cross is defined by the pulse duration r. 
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FIG. 2. The momentum distribution of atoms originally pre- 
pared in state |2) after they interact with a pulse of Raman 
beams, (a), (c) The population transferred to states |1), |3) 
through Raman transitions |2) o |1) and |2) ^ |3), respec- 
tively; (b) cross appearing in state |2) due to these transfers. 
The duration of Raman pulse is r = 2, the Rabi frequency of 
two-photon resonance is a = 0.5. The detunings from ground 
states are = (m = 1, 2, 3). 



As seen from Eqs. (p!Q|) and (p!3|) . the largest excitation 
probability is reached when 

^{Sl-S2)-Pz^Py^'^ = or ^{S3-62)^Pz^Py^l = 0. 

The position of the cross can be changed by adjusting the 
detunings from the ground states, allowing one to excite 
atoms with any two projections of the atomic momentum. 
The excitation of atoms through two Raman transitions 
at once essentially increases the efficiency of cooling in 
a tripod configuration. Note that for Fig. [2] we have 
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used ar = 1 but the short-interaction hmit gives a good 
account for the basic aspects of the process as the relevant 
equations are modified only slightly for increasing ar (see 
Sec.HrCl). 

Similarly to ID Raman cooling, an elementary cycle of 
2D Raman cooling consists of two steps. In the first step, 
atoms prepared in state |2) are selectively transferred by 
a pulse of Raman beams. In order to suppress the Raman 
transfer of atoms with the momentum projection Pz^Py = 
0, we can choose the pulse duration 
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where 6i — S2 = Ss — S2 = S. In the second step, atoms 
return to state |2) via optical pumping. The 7r-polarized 
laser is switched off, the (j_ lasers are tuned into 
resonance with transitions |1) ^ |4) and |3) ^ |4), re- 
spectively. An atom is excited to state |4) with a gain of 
momentum along the z direction, so = pz ^1. Then 
the atom returns to state |2) emitting a photon with mo- 
mentum Ap where \Ap\ = 1. Because of momentum 
conservation, the atomic momentum changes by —Ap. 
So, the atomic population in state |2) becomes 

W2{Pz,Py)f = \a2{Pz,Py)f + \ai{pz-l^Apz,Py^Apy)\^ 

+ |a3(p.+l+Ap',,p^+Ap;)p. 

Different momentum groups can be excited from state 
1 2) by adjusting the duration r and the difference of de- 
tunings S in consistent with (p!5|) . The set of detunings S 
used by us is 



0,1, ..,5. 



(16) 



The position of the cross for large r is situated in im- 
mediate proximity to zero projections, which gives as a 
result that only very cold atoms are left in state |2). 



C. Arbitrary Rabi frequency solution 

To simplify the task, we assume that Si = 63 = S2 -\- 
and consider the only atoms of momentum projection 
Pz = 0. These atoms are characterized by the coherence 
between states |1) and |3), which in turn is derived from 
the difference of Eqs. ([Taj) and ([7c|) 
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Before the Raman pulse starts, the atoms are contained 
in state |2). Hence 61,63 = 0, and one obtains from Eq. 
(fT7|) that 



61 63. 



(18) 



After the substitution of Eq. (p!8|) . the equations ([7]) are 
reduced to a system of two equations 

ibi = (1 + + (5i + 2a) hi + 0^62, (19) 

ih2 = {{Py - if + (^2 + Q^) 62 + 2a6i. (20) 



Then the interaction with Raman beams is described by 
the effective Hamiltonian and probability function 

(21) 

where ^eff = S -\- 2py + a, aeff = ^/2a. 

The Hamiltonian (|2T]) in the dressed-state picture is 
presented by the eigenstates of the probability ampli- 
tudes 



= V^icos^- V^2sin^, ^/4a^^TS^ - Sef^ 

tjan (p — , 

ijj- = ipi sin + V^2 cos ^, 2aeff 

(22) 



with the corresponding eigenfrequencies 



where 7/^2 are the probability amplitudes of states in 
presentation (|2T]) . Expression (|22]) gives original after 
the substitution of ijj^ = 0, = i-^-^ probability 
amplitudes are written as 

= -hi sin (t) e-'^+\ = hi cos e"^^-^. (23) 

By inverting Eq. (|22]) we get 

^1)1 = cos sin 0, iIj2 = —'0+ sin (j) -\- i/j- cos 0, 

which gives via Eq. (f23|) the probability amplitudes 



V^i = y sin2(/) (e" 
V^2 = 62 (sin^0e~ 
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Thus the populations in states |1), |2) in consistent with 
Eqs. (|3]) are given by 



|ai(-l,p^+l)|' = |a0(0,p^)|2 sin'Z^effT, 

^eff 

|a2(0,p^)|2 = \al{{),py)\^ fl - ^ sin^I^effr) , 



eff 



where 20^^ = Y'^(^>2p7+2T~^P^r8c?. Atoms in the 
zero velocity group remain in state |2) for the pulse du- 
ration 
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(24) 



which corresponds to |a2(0,0)p = |a2(0,0)p. One can 
easily see that in the limit of small a this result agrees 
with Eq. ([I5]). 
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FIG. 3. The momentum distribution of atoms after applying 
420 elementary cycles of 2D Raman cooling. Atoms have been 
piled up in a narrow peak near zero momentum with effec- 
tive temperature equal to O.OlTrec- The height of the peak is 
about 80 times larger than that of the original distribution, 
the width (FWHM) is about 0.24 in units of Hk. 



D. Full cooling process 



where /c^ is the Boltzmann constant, M is the Ne atomic 
mass. The duration of only Raman pulses in our scheme 
is about 400 cj^^, where uo^ for metastable Ne equals 
6 /is. Full Raman cooling in turn takes a few millisec- 
onds. The number A/'decay of atoms spontaneously de- 
caying from the upper state to a ground state during a 
Raman pulse is given by 

A^decay = 7^4, A^decay ^ ^^^2 = ^^J^"^^^ 

where is the number of atoms, 7 is the decay rate; 
the number A^4 ^ Ng'^ / /S? of atoms in the upper state is 
obtained from Eq. (|6]). Decay from the upper state can 
be neglected when A/'decay <^ Ni^^. Using Eqs. ([12]) and 
([Hj) one gets the inequality 



iVdecay < 

which can be written in two alternative forms: 



ar > 
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(25) 



(26) 



where in the latter form it was assumed that r ~ 1 and 
a = G'^/uj^A (see Eq. ([5])). Thus by increasing either 
the detuning related to the upper state or the pulse in- 
tensities we can avoid the effects of spontaneous emission 
during the Raman pulses. 



Full Raman cooling consists of varying Raman pulses, 
decreasing the detuning S as in set (p!6|) . After each six 
elementary cycles, the direction of 7r-polarized beam is 
alternated, so the next cycles excite another side of the 
velocity distribution. During cooling process, atoms are 
piled up in a narrow peak near zero momentum. As the 
number of Raman cycles increases, the number of atoms 
in the peak increases, so the momentum distribution be- 
comes narrower, and the temperature of atoms decreases. 
The efficiency of growth increases as the Rabi frequency 
a of two-photon resonance increases. Numerical calcula- 
tions show that deep Raman cooling is only achieved in 
the case of quite large ar. Figure [3] shows the result af- 
ter applying 420 elementary cycles of 2D Raman cooling 
with a = 1. Th e velo city spread of atoms is defined as 
a = (FWHM)/V8 In 2 in units of the recoil velocity Vrec 
It has been reduced from 3vrec to O.l'^rec, corresponding 
to effective temperature Teff = O.OlTrec- 

Consider the 2D Raman cooling under conditions of 
experiment ^] where the same atomic and field con- 
figurations were realized in metastable Ne beam. The 
cr+- and cr_ -polarized waves were provided by lasers with 
^Pump = 588 nm, while laser with Xstokes = 617 nm 
generated the 7r-polarized wave. Hence the final effective 
temperature of cooled atoms is given by 



III. CONCLUSIONS 

We have demonstrated that one can obtain very nar- 
row velocity distributions (<j < 0.1'Urec) simultaneously in 
two directions with our straightforward scheme. The tri- 
pod approach presented here sets a further restriction on 
the available atomic structure, but this is not necessarily 
a major problem. Appropriate optically coupled J = 1 
and J = states are available e.g. in Rb and Ne*. In 
addition to 3D cooling by thermalization, the approach 
can be used to cool atoms in an ID optical lattice, where 
the sample is "sliced" into 2D pancake-like structures. 
This approach is used for investigating optical frequency 
standards with alkaline-earth atoms [6]. The simplicity 
of our approach, coupled with the expected efficiency, 
opens new possibilities for all-optical studies of ultracold 
atoms. 
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